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ABSTRACT 

The non-circular streaming motions in barred galaxies are sensitive to the mass of the bar and can be 
used to lift the degeneracy between disk and dark matter halo encountered when fitting axisymmetric 
rotation curves of disk galaxies. In this paper, we present photometric and kinematic observations of 
NGC 4123, a barred galaxy of modest size {Vrot = 130 km s^^, L = 0.7i»), which reveal strong non- 
circular motions. The bar has straight dust lanes and an inner Lindblad resonance. The disk of NGC 
4123 has no sign of truncation out to 10 scale lengths, and star-forming regions are found well outside 
i?25- A Fabry-Perot Ha velocity field shows velocity jumps of > 100 km at the location of the dust 
lanes within the bar, indicating shocks in the gas flow. VLA observations yield the velocity field of the 
H I disk. Axisymmetric mass models yield good fits to the rotation curve outside the bar regionfor disk 
/-band M/L of 2.25 or less, and dark halos with either isothermal or power-law profiles can fit the data 
well. In a companion paper, we model the full 2-D velocity field, including non-circular motions, to 
determine the stellar M/L and the mass of the dark halo. 

Subject headings: galaxies: kinematics and dynamics — galaxies: halos — galaxies: structure — 
galaxies: individual: NGC 4123 



1. INTRODUCTION 

The rotation speed of gas in a disk galaxy is directly re- 
lated to its mass, providing a dynamical measure of galaxy 
mass if the gas is on circular orbits in the disk plane {e.g. 
Opik 1922; Roberts 1969; TuUy & Fisher 1977). The ro- 
tation curve Vc{R) yields the radial profile of the enclosed 
mass M{< R), modulo the effects of the flattening of the 
mass distribution. Many galaxies have nearly flat rotation 
curves, requiring dark mass more extended than the lumi- 
nous material {e.g. Babcock 1939; Rubin, Ford, & Thon- 
nard 1978; Bosma 1981). 

Unfortunately, even a high quality rotation curve ex- 
tending well beyond the optical disk does not constrain 
the relative contributions of dark and visible matter to 
the centripetal acceleration; such decompositions are al- 
most totally degenerate. For example, van Albada et al. 
(1985) showed that the rotation curve of the well-observed 
galaxy NGC 3198, the "poster child" for dark matter (DM) 
halos, can be fitted equally well with a massive disk and a 
dark halo significant only in the outer parts of the galaxy, 
a less massive disk and a more dominant halo, or even just 
a halo and a zero-mass disk. 

It is critical to try to determine the relative masses of 
disks and halos in galaxies because uncertainty in this 
quantity impacts almost every area of galaxy formation 
and dynamics. It hampers our understanding of: the 
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structure and formation process of disk galaxies; the effi- 
ciency of assembly of baryons into disks; the TuUy-Fisher 
relation between luminosity and rotation width; instabili- 
ties which give rise to spirals and bars, warps, mergers, etc. 
Furthermore, ignorance of dark halo densities and masses 
leaves us without a probe of the initial density fluctuation 
power spectrum on galaxy scales. 

One approach to this problem has been to place a lower 
bound on the mass of the inner DM halo by constructing a 
"maximum disk" model to fit the rotation curve (Kalnajs 
1983; van Albada & Sancisi 1986). Such a model attributes 
the highest possible mass-to-light ratio {M/L) to the lu- 
minous matter, usually consistent with a non-hollow DM 
halo. There are a number of indirect dynamical arguments 
suggesting that the disk mass is in fact close to maximum 
for most bright galaxies {e.g. van Albada & Sancisi 1986; 
Freeman 1992; Debattista & Sellwood 1998) or nearly so 
(Athanassoula, Bosma & Papaioannou 1987; Bosma 1999). 
Counter-arguments against maximum disks have been ad- 
vanced by Efstathiou, Lake, & Negroponte (1982), Bah- 
call & Casertano (1985), Bottema (1993, 1997), van der 
Kruit (1995), Courteau & Rix (1999) and others. Several 
recent studies, motivated by cosmological structure forma- 
tion simulations, have addressed dark halo structure and 
the relation between dark halo size and galaxy size {e.g. 
Navarro et al. 1996, Syer & White 1998, Syer, Mao, & Mo 
1999) and have generally favored strongly sub-maximal 
disk models. 

While models with zero mass disks can fit rotation 
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curves, all reasonable models attribute some mass to the 
luminous component. The disagreement between maxi- 
mum disk and halo-dominated models such as those of 
Bottema (1997) is roughly a factor of 2.5 in disk M/L. 
Population synthesis estimates of stellar M/L ratios can- 
not resolve this discrepancy, because most of the light is 
emitted by massive stars while most of the mass is in faint 
stars (Worthey 1994; Chariot, Worthey & Bressan 1996). 

Barred galaxies offer a way to break the disk-halo de- 
generacy. The elongated bar drives non-circular streaming 
motions in stars and gas within the bar region. The mag- 
nitude of the non-circular motions depends mostly on the 
mass of the bar and, to a lesser extent, on the pattern 
speed at which it rotates. The departures from axisym- 
metry in both light distribution and kinematics can be 
observed. The (unobservable) dark halo is thought to be 
mostly "pressure" supported and therefore more rounded 
than the bar. The stellar bar is therefore the dominant 
source of non-circular streaming motions, which allow us 
to estimate its M/L, and by extension that of the stellar 
disk. 

In practice, we need extensive observations and model- 
ing to carry out this program. The most easily observable 
kinematic tracers are optical and radio emission lines from 
gas in the galaxy disk. Because the gas is not on circular 
orbits, and we observe only line-of-sight velocities, we can- 
not infer the full space velocity field of the gas, as we can 
for an axisymmetric galaxy. Furthermore, gas flow in a 
barred potential generally gives rise to shocks, which pre- 
vent us from directly inverting the observed gas velocities 
to derive the mass distribution. We therefore need to de- 
rive the bar mass by matching fluid models incorporating 
shocks to the observed velocity field. 

In this paper we describe a detailed study of the barred 
galaxy NGC 4123. We selected this galaxy for its fairly 
strong bar, a moderate inclination, favorable projection 
angle of the bar, colors and morphology suggesting sub- 
stantial Ha emission, and modest size. Because projection 
of a non-axisymmetric velocity field into the observable 
line-of-sight velocity depends on its angle to the galaxy 
major axis (line of nodes), candidates for study must be 
selected carefully (Pence & Blackman 1984b; Long 1991). 
If the bar is aligned with either the major or minor axis 
of the galaxy disk as it appears on the sky {e.g. NGC 
1300, England 1989), it is difficult to separate the circular 
and non-circular parts of the velocity field after projection, 
since the orbits or streamlines are elongated in approxi- 
mately the same direction as the bar {e.g. Athanassoula 
1992b). The projection is most favorable for bars at inter- 
mediate position angles - a 45° angle between the bar and 
the line of nodes is best. 

We require two-dimensional surface photometry in sev- 
eral broad-band filters and a 2-D velocity map over the en- 
tire galaxy. "Three-dimensional" imaging spectroscopy is 
possible with Fabry-Perot interferometers at optical wave- 
lengths, and with aperture-synthesis interferometers at ra- 
dio wavelengths. These techniques are complementary in 
angular resolution and field of view, and the spectral lines 
arise from different phases of the interstellar medium. This 
galaxy is the first of several barred galaxies we have ob- 
served with the Rutgers Imaging Fabry-Perot Interferom- 
eter at CTIO (Schommer et al. 1993). 

We present surface photometry and velocity maps of 



Parameter 



Value 



Recession velocity 
Assumed distance 
ms (within 3i?25) 
my (within 3i?25) 
mi (within 3i?25) 
PA 

inclination 

Li (extinction corrected) 

Vrot 

21 cm fiux 
M„i 



1327 km 

22.4 Mpc 
12.06 
11.45 
10.33 
-33 ±2° 
45 ±4 ° 
1.72 X 10^° L 
130 km s-i 

63.5 Jy km s~^ 
7.5 X 10^ Mq 







Table 1 
NGC 4123 - Basic Properties 



NGC 4123, from which we construct axisymmetric models. 
Mass models with a maximum disk M/L of 2.25 and either 
isothermal or power-law profiles can fit the rotation curve 
well. In a companion paper (Weiner et al. 2000, hereafter 
Paper II) , we model the full 2-D velocity field to determine 
the separate masses of the disk and halo. 

2. NGC 4123 - BASIC PROPERTIES 

NGC 4123 is of moderate size and has no pronounced 
morphological peculiarities. It is an SB(r)c galaxy with 
total blue magnitude Bt = 11.98, a rotation amplitude 
of approximately 130 km s^^, and a recession velocity of 
1327 km s~^, from the RC3 catalog (de Vaucouleurs et 
al. 1991), via NED0. We adopt a distance of 22.4 Mpc 
(m — M = 31.75), based on a Hubble constant of 75 
km Mpc~^ and a correction to the microwave back- 
ground rest frame, with no Virgocentric correction. Dy- 
namical estimates of M/L vary inversely with the adopted 
distance D; M/L values given in this paper are in so- 
lar units and are implicitly followed by the factor /175 
{Ho/75 km s~^ Mpc~^). However, estimates of the de- 
gree of disk maximality (disk-to-halo ratio) do not depend 
on adopted distance. The basic properties of NGC 4123 
are summarized in Table ^ 

NGC 4123 has a fairly strong bar, but the bar does not 
totally dominate the disk (unhke e.g. NGC 1097 or NGC 
1365; cf. Ondrechen & van der Hulst 1989; Ondrechen 
et al. 1989; Lindblad, Lindblad, & Athanassoula 1996), 
so it is not an extreme case. The galaxy has a favorable 
inclination (45° as determined from 21 cm observations, 
see below) and the bar is aligned at approximately 53° to 
the line of nodes when deprojected, which is almost ideal 
for recovering non-circular motions. 

NGC 4123 has a companion galaxy, NGC 4116, which 
lies 14' southwest of NGC 4123 and at nearly the same 
systemic velocity, 1312 km s~^. NGC 4116 is of SB(rs)dm 
type and is fainter than NGC 4123, with magnitude Bt = 
12.41, from the RC3 via NED. At 22.4 Mpc, the separation 
of the galaxies in the plane of the sky is just 91 kpc. How- 
ever, neither shows much sign of tidal interaction (see Sec- 

^The NASA/IPAC Extragalactic Database (NED) is operated by 
the Jet Propulsion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and Space Adminis- 
tration. 
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tion 1^) and it is possible that the real separation is much 
greater. The largely undisturbed kinematics of the two 
galaxies, shown in Taylor et al. (1995) and discussed fur- 
ther below, indicate that any interaction does not visibly 
affect their internal dynamics. The conclusions we draw 
from studying the dynamics of the bar region of NGC 4123 
are unaffected by the presence of NGC 4116. 

3. SURFACE PHOTOMETRY 

We require surface photometry in order to estimate the 
distribution of luminous mass in the galaxy. It is usual to 
assume that mass is proportional to luminosity; a uniform 
M/L has the advantage of being the simplest possible as- 
sumption, and is generally thought to be reasonable when 
color variations are small. This is difficult to prove, since 
the bulk of the mass is in low mass stars while high mass 
stars emit much more light per unit mass. Furthermore, 
extinction is particularly variable in barred galaxies, which 
generally have strong dust lanes within the bar. Such 
problems are lessened at red and near IR wavelengths, 
since bluer optical bands are more strongly affected by 
emission from young massive stars and by spatially vari- 
able extinction, although / (and even K) band observa- 
tions are still affected by young stars and by dust (de Jong 
1996b, Rhoads 1998). 

3.1. Observations 

The initial imaging observations of NGC 4123 were 
made on the night of 1994 April 1 with the 0.9-meter tele- 
scope of Cerro Tololo Interamerican Observatory.^ We 
used a Tektronix IK x IK CCD at the f/13.5 Cassegrain 
focus, with a pixel scale of 0.39", and a field of view of 
6.6'. We observed NGC 4123 in BVRI filters, obtaining 
10 minute exposures in BVR and 3 x 10 minute exposures 
in /, offset between exposures to improve flatness of the 
final image. The seeing was 1.2" in VRI, or 130 pc at the 
adopted distance of 22.4 Mpc, and 1.4" in B. We reduced 
the images with IRAF^, subtracting overscan and a bias 
frame made by a robust combination of many individual 
biases, and flatfielding with twilight flats in the usual way. 
Photometric standard stars from the E-regions (Graham 
1982) were observed at several times during the night and 
used to derive extinction coefhcients. 

The field of view of the CTIO observations is rather 
small for such a large, nearby galaxy. We reobserved NGC 
4123 to obtain a larger field, using the Swope 1-meter tele- 
scope at Las Campanas Observatory on the night of 1998 
January 25. These observations used the LCO Tek#5 2K 
X 2K CCD at the f/7 Cassegrain focus, with a pixel scale 
of 0.70" and field of view of 24'. We obtained 10 minute 
exposures of NGC 4123 in BV filters and 2 x 10 minute off- 
set exposures in the / filter. The seeing was approximately 
1.6" in VI and 1.8" in B. The images were again reduced 
with IRAF, subtracting overscan and bias and flatfielding 
with twilight flats. Photometric standard stars from Lan- 
dolt (1992) were observed throughout the night and used 
to derive extinction coefficients. The photometric calibra- 

^CTIO is operated by the Association of Universities for Research 
in Astronomy, Inc. under a cooperative agreement with the National 
Science Foundation. 

^IRAF is distributed by NOAO, which is operated by AURA Inc., 
under a cooperative agreement with the National Science Founda- 
tion. 



tions of the CTIO and LCO data are consistent. 

The LCO Tek#5 chip exhibits fringing and illumination 
gradients in the / filter which cannot be removed with 
twilight flats, since the night sky spectrum is of a different 
color and contains strong emission lines. We constructed 
a supersky image by combining a number of dithered / 
images from throughout the 3-night run, scaled by the sky 
brightness, and subtracted this image from the NGC 4123 
/ image. This process removed both small-scale fringes 
and large-scale gradients in the image. 

3.2. Morphology 

The inner part of the combined CTIO 0.9-m / image 
of NGC 4123 is presented in Figure |l|. The bar has a 
semi-major axis of approximately 50" on the sky, as de- 
termined by the transition from boxy to essentially ellip- 
tical isophotes, and by the break in the surface brightness 
profile (see below), or 5.6 kpc when deprojected. The bar 
has a rather complex shape, with a bright boxy compo- 
nent near the center and a fainter elongated component, 
which is of constant width or even broadens a little to- 
wards the ends; there is some evidence for an isophotal 
twist between the two components. The complex shape 
makes it difficult to define a single axis ratio for the bar: 
the elongated component has an axis ratio of about 6:1, 
while the boxy component has an axis ratio of about 2.5:1. 

Other than the central source (discussed below), NGC 
4123 has no significant spherical bulge, which makes it 
particularly suitable for the deprojection method we use in 
our simulations of the gas flow in this galaxy (see Paper II) . 
Outside the bar region, the disk of NGC 4123 is relatively 
smooth in / light, with some spiral structure but no sharp 
features or strong brightness variations. 

There is an isophotal twist in the outer regions of the 
disk. At very large radii the PA and axis ratio of the 
isophotes are consistent with the PA of the galaxy major 
axis and the inclination as derived from 21 cm observations 
of the H I disk of the galaxy (see section ||). However, at 
radii less than ~ 100", but well outside the bar, the PA 
of the disk isophotes is intermediate between the PA of 
the bar and the PA of the kinematic major axis indicated 
by the 21 cm observations. A heavily smoothed image 
shows a faint spiral arm pattern in the outer disk, which is 
responsible for the isophotal twist. There is no large-scale 
twist in the kinematic line of nodes in this range of radii 
in the H I observations, which indicates that the isophotal 
twist is not accompanied by a warp (Briggs 1990). 

The B — I color map, presented in Figure |2| was con- 
structed from CTIO B and / images convolved to 1.8" res- 
olution. The dust lanes, conspicuous in this color map, are 
linear features, offset from the bar major axis towards its 
leading edge. Dust lanes of this geometry are common in 
strongly barred galaxies (Sandage 1961; Sandage & Bedke 
1988; Athanassoula 1992b and references therein). Color- 
color plots of pixels within the bar (as in Quillen et al. 
1995) show that in the east dust lane, the slopes of i? — / 
and V — I versus B — I are consistent with a dust absorbing 
screen with Ry ^ 3.1 (Cardelh, Clayton & Mathis 1989), 
confirming that the dust lane in the east half is near to us. 
In the west dust lane the slopes are shallower, so that the 
dust lane is on the far side of the bar, as expected. 

The small bright nucleus at the center of the galaxy 
stands out as the central blue spot. It is barely resolved 
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Fig. 1.— NGC 4123, / image. 

The inner 240" square region of tlie CTIO 0.9-m /-band image of NGC 4123, on a linear scale. North is up and west is to the right. The major axis of 
the disk is 33° west of north. [Resolution degraded by a factor of 2 for preprint version.] 



by these observations. Its color, B — I = 1.4, is quite 
blue compared to its surroundings [B — I = 2.2), and 
its luminosity in the / band is 1.8 x lO^i©. The center of 
NGC 4123 appears in lists of known emission line and UV- 
excess sources as UM 477 (MacAlpine & Williams 1981) 
and Markarian 1466, suggesting that it is a nuclear star- 
burst. It is likely, therefore, that the M/L of the nucleus 
is lower than the overall M/L of the stellar disk, a point 
we return to in the models of Paper II. 

3.3. Surface brightness profiles 

We constructed axisymmetrized surface brightness pro- 
files of NGC 4123 from the Las Campanas BVI data. The 
21' field allows us to extend the surface brightness profiles 
to large galactic radius. 

We first masked approximately 1100 stars and faint 
galaxies and a small number of cosmic rays out of the 
image, using large masks around the stars to exclude the 



faint halo of the stellar PSF. The surface brightness pro- 
files are then constructed by dividing the image into el- 
liptical annuli, 1.7" wide along the major axis, with the 
disk inclination of 45° and position angle of -33° as de- 
rived from the 21cm H I observations (see below). In each 
annulus, we compute the trimmed mean with a 2-sigma 
clipping. We determined the sky level by computing the 
biweight, a robust estimate of the mode (Beers, Flynn & 
Gebhardt 1990), of all pixels within an elliptical annulus 
of the same inclination and PA, an inner major axis radius 
of 551" (60 kpc) and outer major axis radius of 593" (65 
kpc). The error estimate is computed directly from the 
dispersion of the pixel values in each annulus, excluding 
clipped pixels, with the statistical error of the sky value 
added in quadrature. 

The axisymmetrized i?, V , and / surface brightness pro- 
files of NGC 4123, and the color profiles, are shown in 
Figure 0. We have not included a correction for internal 
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Fig. 2.— NGC 4123, B - I color map 
The inner 120 square region of the CTIO B — I color map of NGC 4123. North is up and west is to the right. Darker colors arc redder. The dust 
lanes in the bar arc apparent as the dark lines running east and west of the nucleus. The nucleus of the galaxy is much bluer than its surroundings. 
The blue spots at the west end of the bar are foreground stars. [Resolution degraded by a factor of 2 for preprint version.] 



extinction in these plots, but we have corrected for geo- 
metric inchnation to obtain a face-on surface brightness: 

Mcorr = Mobs + 2.5 log(b/a), (1) 

as appropriate for a transparent galaxy. At an inclination 
of 45° this correction is -1-0.38 mag arcsec"^. The shape of 
the profiles, the colors, and the color gradient are typical 
of SBbc galaxies. 

There is a substantial color gradient, ~ 0.4 mag, inside 
the bar radius, but this does not imply an equally large 
color gradient in the old stellar population. Inside the bar, 
the strong dust lanes produce reddening, while outside the 
bar, many H ii regions with very blue colors are present 
in the spiral arms. We measured mean colors in several 
regions in the bar and disk in and out of dust lanes and 
H II regions. The results are given in Table |. When 
obvious dust lanes and H ii regions are excluded, there is 



Region 


B — I (mag) 


Dust lane in bar 


2.3 


Bar/lens outside dust lane 


2.0 


Disk, interarm region 


1.8 


Disk, spiral arms 


1.5 



Table 2 
NGC 4123 - Mean Colors 



only an 0.2 mag difference in B — I between the bar and 
disk, so there is no evidence for a significant difference in 
stellar population. Since unresolved dust lanes and knots 
of young stars may be present, the actual color difference 
could be smaller. 

In our subsequent use of the / data to fit the rotation 
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curve and model the potential, we apply a —0.15 mag over- 
all correction for internal extinction in the / band, calcu- 
lated from the Amt — — 1.01og(&/a) prescription given by 
Giovanelli et al. (1994), which Palunas (1996) found to be 
reliable for similar /-band data. Galactic extinction Acxt 
towards NGC 4123 is negligible (Burstein & Heiles 1982). 

We did not attempt to correct the / image for differen- 
tial extinction within the galaxy. This procedure requires 
extensive multicolor data, including near-IR imaging, and 
detailed modeling of radiative transfer (see Byun et al. 
1994; Regan et al. 1995; de Jong 1996b). To estimate the 
extinction in the east (stronger) dust lane, we selected a 
region of the bar where the dust lane appeared strong in 
B and took a cut through the data perpendicular to the 
dust lane. Assuming the unextincted bar light profile is 
symmetric when reflected about the bar major axis, we 
estimate the differential extinction to be ~ 25% (0.3 mag) 
in B and ~ 10% (0.1 mag) in /. Alternatively, for the dust 
lane reddening of Eb-i = 0.3, an Ry = 3.1 extinction law 
(Cardelli, Clayton & Mathis 1989), and an absorbing dust 
screen, the dust lane extinction in / is 0.17 mag. Since the 
dust lanes cover only a small part of the bar, the amount 
of / light lost to extinction is just a few percent. 
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Fig. 3. — BVI surface brightness and color profiles of NGC 4123 
Upper panel: The axisymmctrizcd LCO B, V, 7-band surface brightness 
profiles of NGC 4123, in magnitudes arcsec""^. The profiles have been 
corrected for inclination but not internal extinction. B, I are triangles, 
squares, and circles respectively. The /-band profile is plotted to 40 kpc 
while B and V arc plotted to 35 kpc. 

Lower panel: The axisymmctrizcd B — V and B ~ 1 color profiles of NGC 
4123, in magnitudes arcscc^^. The profiles have not been adjusted for 
internal reddening. 

The luminosity profile has a strong peak at the center, 
a shoulder in the bar region, and an approximately expo- 
nential decrease outside the bar, though the slope changes 
somewhat. The radius at which the profile breaks to ex- 
ponential, ^ 50", or 5.5 kpc, is the same as the bar length. 
This type of surface brightness profile, with flat and expo- 
nential regions, was designated type II by Freeman (1970), 



and is fairly common in both axisymmetric and barred 
galaxies; see for example the sample of Palunas & Williams 
(2000). 

Table |^ summarizes overall photometric parameters of 
NGC 4123 from the LCO observations. The nominal opti- 
cal radius of the galaxy (i?25 in B) is 11.1 kpc; R23.5 in I is 
11.3 kpc. The exponential disk scale lengths given in Table 
m were determined from least-squares fits between 6 and 
10 kpc, where the data are well described by an exponen- 
tial. The inward extrapolated central surface brightness in 
/ of 19.59 mag (19.44 mag when corrected for extinction) 
makes this galaxy relatively high surface brightness, but 
not extremely so (de Jong 1996a; Palunas 1996). The total 
/ magnitude within 3i?25 is 10.33; the absolute magnitude, 
corrected for internal extinction, is Mi — —21.57 and lu- 
minosity L/ = 1.72 X 10^° Lq. The absolute B magnitude, 
uncorrected for extinction, is Mb = —19.7, a luminosity 
of 0.7i,. 

3.4. Outer disk profiles 

We have been able to trace the surface brightness pro- 
files to quite large radii, beyond 2i?25, even to 3i?25 for 
/ (9 magnitudes below the sky). At very large radii the 
profiles become uncertain due to a combination of photon 
noise and errors in the sky level. Beyond 20 kpc, the color 
profiles are probably not reliable - in particular, the red- 
dening trend seen in B — I but not in B — V could be due 
to a systematic error in the / sky level. Some small-scale 
features in the profiles are certainly real, such as the red 
bump in B — / at 12-14 kpc and the bump in all three fil- 
ters at 16 kpc, which both correspond to faint spiral arms. 
There is also a spiral arm at ~ 24 kpc. 

The surface brightness profiles show no sign of reaching a 
plateau before they disappear into the noise, at i? 30 kpc 
for B and V, and 35 kpc for /, so the large extent of 
the galaxy is not due to undersubtraction of the sky. We 
constructed an extended PSF from observations of bright 
stars to measure the contribution of scattered light, which 
is negligible. Scattered light outside a sharply truncated 
disk would produce a surface brightness profile with an 
observed drop of 4 magnitudes at the disk edge, which is 
clearly ruled out. 

In all three filters, the galaxy is easily visible to i? 
25 kpc when displayed at high contrast, and to greater 
extent when the images are smoothed with a median filter. 
Spiral arms are traceable to i? ~ 25 kpc; the arms are 
strongest in the B filter, while in the / filter the extended 
galaxy light is more smoothly distributed than in B and 
V. 

There is no evidence in the data that the disk is trun- 
cated at large radius, in sharp contrast with results that 
suggest disk galaxies are generally truncated at 3 to 5 scale 
lengths {e.g. van der Kruit & Searle 1981a, b, 1982; van der 
Kruit 1989; Barteldrees & Dettmar 1994; Morrison, Boro- 
son & Harding 1994). The SB profiles of NGC 4123 are 
traceable twice as far, to 10 scale lengths in the / filter. 
At least 20% of the galaxy luminosity is outside R25. One 
other galaxy, NGC 5383, definitely does not have a trun- 
cated disk (Barton & Thompson 1997). 

The origin of this extended stellar disk is uncertain. It 
has been suggested that truncated disks occur because a 
stability condition cuts off star formation below a critical 
gas surface density (Kennicutt 1989; Sellwood & Balbus 
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V/T /t\1 "Vi^l 












brightness 




Mo 


scalelength 


R25 


2i?25 3i?25 




mag/n" 


mag/n" 


mag/n" 


kpc 




mag 


B 


22.35 


25.11 


21.38 


3.46 


12.29 


12.09 12.06 


V 


21.52 


24.64 


20.65 


3.25 


11.65 


11.46 11.45 


I 


19.43 


23.38 


19.59 


3.20 


10.55 


10.36 10.33 



Table 3 

NGC 4123 - LCO Photometric Data 



Note. — The central surface brightness /ig in column 4 has been 
corrected for inclination, but not for internal extinction. 



1999), and it is possible that stars from the inner disk 
could be scattered to large radius by a bar-like perturba- 
tion (Schwarz 1984). However, NGC 4123 has faint spiral 
arms extending to ~ 25 kpc, with small blue knots that are 
likely sites of massive star formation. These star-forming 
regions are similar to the H II regions in faint spiral arms 
that Ferguson et al. (1998) detected to 2i?25 and beyond, 
in deep Ha imaging of three late- type spirals. Thus at 
least some of the extended disk of NGC 4123 formed at 
large radius, although massive star formation in the outer 
disk appears to be confined to spiral arms, as discussed by 
Ferguson et al. (1998). The blue knots have luminosities 
from Mb = —12.5 to at least 2 magnitudes fainter, and 
mean colors B — V ^ 0.1, V — I ^ 0.0. The very blue colors 
indicate that the knots are younger than 10^ yr (Bruzual 
& Chariot 1993). The largest knot has Lb = 10^ Lq, sug- 
gesting M ~ 2 X 10^ M0. The star formation rate implied 
outside R25 is consistent with the range found by Ferguson 
et al. (1998). 

4. FABRY-PEROT IMAGING SPECTROSCOPY 
4.1. Observations 

We observed NGC 4123 in Ha emission with the Rutgers 
Imaging Fabry-Perot Interferometer to obtain a 2~D ve- 
locity field with high spatial resolution. The observations 
were made at the Cassegrain focus of the CTIO Blanco 
4-meter telescope, on the night of 1994 April 3. We used a 
Tektronix IK x IK CCD detector, with 0.36" pixels and a 
field of view of 160". The Fabry-Perot etalon has spectral 
resolution of 2.5 A at Ha, yielding a velocity resolution of 
FWHM ~ 150 km s~^. The instrumental spectral profile 
is very well approximated by a Voigt function, with Gaus- 
sian ac = 26 km s~^and Lorentzian cr^ — 33 km s~^. The 
free spectral range of this etalon is 85 A at this wavelength, 
and a blocking filter of 75 A FWHM ensured that only one 
order was transmitted. 

The Fabry-Perot instrument provides a normal spa- 
tial image of the sky, but with a very narrow, tunable, 
bandpass. The image is not monochromatic, but varies 
quadratically with radius from the optical axis of the in- 
strument near the center of the image, with a 6 A gradi- 
ent. There are two temporal variations in the instrument: 
the wavelength zeropoint of the etalon varies due to tem- 
poral drift of the control electronics, and the position of 
the optical axis changes due to spectrograph flexure. In 
practice, these complications cause little difficulty; we take 



hourly calibration exposures while observing to monitor 
the drifts, and the data analysis procedures account for 
them and for the radial wavelength gradient. 

The observations of NGC 4123 are 18 10-minute expo- 
sures spanning a range of 20 A in steps of approximately 
1.2 A (54 km s-i). The images were processed using 
IRAF, and each flatfielded with a dome flat taken at the 
appropriate wavelength. The images were then cleaned of 
cosmic ray events, shifted into spatial register, convolved 
to a common PSF of 1.5" FWHM, and sky subtracted 
to construct a data cube. At each pixel in the image we 
smoothed the data cube spatially with a small Gaussian 
kernel of 0.6" FWHM, and then fitted a Voigt function 
along the spectral axis. The fits yielded maps of the veloc- 
ity, line strength, continuum level, intrinsic line width, and 
estimates of the uncertainties in these quantities. Points 
with intensities below a relatively small threshold or with 
velocity uncertainties greater than 12 km s^^ are deleted 
from the maps. 

The velocity map of NGC 4123 produced from the 
Fabry-Perot observations is presented in Figure ^. The 
orientation is the same as that of Figure |l|, although the 
scale is different. The colormap indicates heliocentric ve- 
locity; the approaching side of the galaxy is blue, and the 
receding side is red. The diagonal tick marks indicate the 
minor axis of the disk. Assuming that the spiral arms 
are trailing, the northeast side of the galaxy is the side 
closer to us. The 160" field of view does not extend much 
outside the bar region. The velocity map from 21-cm ob- 
servations, presented in section ^, extends to much greater 
galactic radius. 

4.2. The streaming motions in the bar 

The non-axisymmetric streaming motions associated 
with the bar can be seen in Figure ^ In an axisym- 
metric galaxy, the contour of the systemic velocity would 
run straight down the minor axis, dividing the map into 
approaching and receding sides. In NGC 4123, the gas 
streamlines elongated along the bar cause the isovelocity 
contours to twist and run more closely along the bar, pro- 
ducing the S-shaped isovels characteristic of barred galax- 
ies (Bosma 1981). To the east and northeast of the galaxy 
center, there are receding velocities (yellow in the colormap 
of Figure ^ in the bar region that are on the approaching 
side of the minor axis, and vice versa to the west. 

Furthermore, the isovelocity contours are not straight 
along the bar, but are offset from the bar major axis. The 
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Fig. 4. — The radial velocity field of the bar region of NGC 4123 
Velocity map of the inner regions of NGC 4123, from Faiiry-Pcrot observations of the Ha emission line. North is up and west is to the right. The 
area shown is 170" diameter, smaller than that of Figure ll Velocity contours arpiplotted over the colormap at intervals of 25 km . The vertical 
line to the east of the center indicates the cut through the aata plotted in Figure pi and in Paper II; the data between the two ellipses are used in the 
goodness-of-fit tests in Paper II. The diagonal lines to NE and SW indicate the minor axis of the disk. 



isovels twist as they pass through the center of the bar. 
For example, the contour of systemic velocity, blue-green 
in the color map of Figure ^, is offset toward north on the 
east side of the bar (left of center in Figure makes a 
jog as it passes through the center, and is offset toward the 
south on the west side of the bar (right of center in Figure 
I). The regions of strong Ha emission in the bar are in 
similarly offset straight lines. There is also a sharp jump 
in velocity across the bar, i.e. parallel to the minor axis of 
the bar. These features suggest the presence of shocks in 
the gas flow pattern, at the location of the dust lanes in 
the bar. 

Dust lanes are presumably caused by a build-up of gas 
at shocks in the gas flow, as first suggested by Prcndergast 
(1966) and by many subsequent authors {e.g. van Albada 
& Sanders 1982; Prcndergast 1983; Athanassoula 1992b; 
but cf. Beck et al. 1999). These shocks are thought to 



be caused by the non-circular gas motions along the bar: 
gas falling down along the bar potential reaches high ve- 
locities, and as the gas climbs away from the center, up 
the potential, it decelerates. Eventually a pile-up occurs 
and the shock is formed. Velocity jumps associated with 
bar dust lanes have been observed, in e.g. NGC 5383 (Du- 
val & Athanassoula 1983), NGC 6221 (Pence & Blackman 
1984a), NGC 1365 (Lindblad et al. 1996a), and NGC 1530 
(Regan, Vogel & Teuben 1997). 

The regions of strong Ha emission aligned with the dust 
lanes in NGC 4123 suggest that compression of the gas by 
the shocks triggers star formation. The offset of the dust 
lanes and the associated jog of the velocity contours as 
they pass through the galaxy center are evidence for an in- 
ner Lindblad resonance (ILR) at the bar center (Athanas- 
soula 1992a,b). An ILR can be inferred when some stel- 
lar orbits near the center are elongated perpendicular to 
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the bar, rather than along it; the gas streamhnes twist to 
follow these stellar orbits {e.g. Figure 2 of Athanassoula 
1992b). This twist creates the jog in the velocity contours 
and the offset in the shocks along the bar. 

The sharp jump in velocity across the bar of NGC 4123, 
plotted in the upper panel of Figure ^ is strong kinematic 
evidence that there are indeed shocks in the bar. The 
radial velocity in the H ii regions about 20" east of the 
galaxy center changes by 80 km s""'^ in just 4". The jump 
is 115 km s~^ when corrected for the galaxy inclination 

- and this is just in the line-of-sight component of the 
velocity; the absolute velocity jump must be larger. There 
are also large velocity gradients in the galactic nucleus, 
inside the ILR, typically 40-50 km s^^ in observed velocity 
(55-70 km s~^ after inclination correction) across just 2". 
This angular size is barely resolved by these observations 
and the actual gradient is likely to be even steeper. The 
velocity gradient shown in Figure || is at the same location 
as the dust lane, confirming that the shock and dust lane 
are associated. 

Because the shocks are closely linked to the non-circular 
gas motions, their properties are dependent on the mass of 
the bar, which controls the ellipticity of the potential. For 
this reason, the strength and location of the shocks will 
prove crucial for measuring the M/L and pattern speed, 
discussed in Paper II. 

There is Ha emission at distances from the galaxy's cen- 
ter of less than the bar radius, but not within the bar itself 

- for example, the small H ii regions some 20" northwest of 
the galaxy center, and the large H ii region 20" southeast 
of the center. These regions have large velocities relative 
to systemic, shown in Figure ^ and provide useful data 
points for constraining the models of Paper II. 

Outside the bar radius, the non-axisymmetric motions 
fall off rather quickly and the contour of systemic velocity 
becomes more closely aligned with the minor axis. There 
are spiral arms with many H ii regions visible to the north- 
west and southeast, but the velocity perturbations associ- 
ated with them are much smaller than those in the bar. 
The limited coverage and field of view of the Ha veloc- 
ity map makes the global axisymmetric velocity field more 
difficult to pick out; the 21 cm H I velocity map, presented 
below, shows it well. 

5. RADIO INTERFEROMETRY 

We have used the VLA to image the neutral hydrogen 
in NGC 4123. The H I data are complementary to the Ha 
data. The higher angular resolution of the Fabry-Perot 
data can trace the kinematics within the bar, while the 
VLA data can image diffuse H I and trace the rotation 
curve far beyond the nominal optical radius i?25- 

5.1. Observations 

We observed NGC 4123 in the H I 21 cm fine with the 
3km (C) and 1km (D) array of the Very Large Array 
(VLA) of the National Radio Astronomy ObservatoryPl 
for 8 and 2 hours in December 1994 and April 1995 re- 
spectively. This combination gives an angular resolution 
of ~ 15", and provided 35 independent beams along the 
galaxy major axis, with sensitivity to smooth structures 

® The NRAO is operated by Associated Universities, Inc., under 
a cooperative agreement with the National Science Foundation. 
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Fig. 5. — Radial velocity and surface brightness profiles crossing 
the bar 

The top panel shows line-of-stight velocity for the north-south line across 
the bar indicated in Figure Ml. South is at negative x\ the direction of 
rotation of the galaxy is fromTeft to right. The bottom panel shows V and 
/ surface brightness profiles along the same line. The surface brightness 
profiles have been sky-subtracted and smoothed with a 3 pixel (l.l") 
boxcar filter. The units are linear (counts/pixel) and the V profile has 
been scaled up by a factor of 3. The position of the prominent dustlane 
is indicated. 



up to 12'. We used 63 channels with a channel width of 
10.5 km s^^. Standard AIPS calibration procedures were 
used and the UV data of the two arrays combined. We sub- 
tracted the continuum in the u — v domain using UVLIN, 
using 19 line-free channels. 

We imaged the data with the AIPS task IMAGR using 
robust weighting (Briggs 1995), yielding a sensitivity al- 
most equal to natural weighting and a resolution almost 
equal to uniform weighting. An elliptical Gaussian fit to 
the beam yielded a beam shape of 19.1" x 16.1" FWHM 
with a position angle of 56.5°. The rms noise in line-free 
channels after continuum subtraction is 0.38 mJy/beam. 
We cleaned the images using the CLEAN algorithm to a 
flux limit of 0.19 mJy/beam, and made images of total H I 
intensity, the intensity weighted velocity field, and the ve- 
locity dispersion with the AIPS task MOMNT. The total 
21 cm fiux observed is 63.5 Jy km s~^, which recovers all 
the single-dish fiux as tabulated in Huchtmcicr & Richter 
(1989). 

5.2. Discussion 

The VLA 21 cm intensity and velocity field shown in the 
left panel of Figure ^ reveals that NGC 4123 has a large H I 
disk of 11' X 7.5', and a regular velocity field. The high- 
est H I column densities occur in a ring where the surface 
brightness is mostly above 10^^ cm~^. The ring's radius is 
some 70", or 7.6 kpc, placing it somewhat outside the bar. 
A loose spiral pattern is also visible outside the ring. In- 
side the ring, the surface brightness decreases toward the 
center, from 5 x lO^'' cm^^ to a low of 4 x 10^^ cm~^ at 
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Fig. 6. — 21 cm emission, velocity field, and residual velocity 
Left panel: The 21 cm surface brightness is plotted as grayscale, with velocity contours overlaid. North is up and west is to the right. The surface 
brightness grayscale runs from Nh — to 10^"^ cm~^. The contours are of heliocentric velocity from 1225 to 1425 km s~^ at intervals of 25 km s"^, 
with receding side in the northwest. 

Right panel: The residuals from our rotation curve extraction arc plotted (data - model), on a grayscale from -20 to 20 km s^^. Non-circular motions 
are seen in the bar and in the outer spiral pattern. 



the center. The hole in the center of the gas distribution 
is frequently seen in barred galaxies; e.g. NGC 1300 (Eng- 
land 1989) and NGC 1398 (Moore & Gottesman 1995). 
It is presumably caused by the bar, which sweeps gas in- 
side its corotation radius from the bar region towards the 
galaxy center, as seen in fluid-dynamical simulations {e.g. 
Athanassoula 1992b; Finer et al. 1995; Weiner & Sellwood 
1999; Paper II). 

The overall velocity field in Figure ^ is quite regular. 
The effect of the bar is visible only inside a radius of ^ 30"; 
the contours close to the systemic velocity in this region 
are twisted to run more east-west, along the bar, rather 
than along the minor axis of the galaxy. 

In the outer parts of the H I disk, velocity perturbations 
due to the spiral arms are visible as kinks in the velocity 
contours. These are particularly noticeable in the velocity 
contours about 160" northeast of the galaxy center. The 
contours bow to the northwest as they cross the spiral arm. 
The effect is of order 10-20 km s^^ in line-of-sight velocity. 
Similar spiral-arm-related kinks of lesser magnitude are 
visible in other regions of the velocity field. 

The only other unusual features in the H I distribution 
are the low-column-density arm or extension on the north- 
west side of the disk, and the even lower column-density 
fragments of H I emission at the north of the disk. The 
shape of these features suggests that they could be due to 
some kind of tidal interaction, although because they are 
close to the detection limit, the "fragments" may not be 
as isolated as they appear. In any case, the kinematics 
of the arm and fragments are completely consistent with 
unperturbed circular rotation. There is no obvious com- 
panion very close to the arm, either in the Palomar Sky 
Survey, the Digitized Sky Survey, our BVI imaging, or the 



21 cm observations. As mentioned in the introduction, the 
galaxy NGC 4116 is just 14' to the southwest of NGC 4123 
and is also visible in our VLA observations. We included 
it in the imaging and CLEANing process and mapped its 
H I emission and velocity. Figure |^ shows the total H I 
image and velocity field of both galaxies, corrected for the 
primary beam response. NGC 4116 has an H I disk just 
over half the angular size of the H I disk of NGC 4123, and 
kinematics consistent with circular rotation and a slowly 
rising rotation curve. 

Tides from NGC 4116 could be responsible for the arm 
of NGC 4123, but if this is the explanation, it is curious 
that the counter-arm pointing away from the companion is 
the more pronounced; furthermore, the larger NGC 4123 
does not appear to have distorted the smaller galaxy. The 
kinematics of the arm can be modeled in circular rotation 
with the rest of the H I disk of NGC 4123, as described 
below. Moreover, there is no sign that NGC 4116 has 
any significant influence on the internal dynamics of NGC 
4123. 

5.3. The rotation curve of NGC 4123 

We derived a rotation curve from the H I data for NGC 
4123 using a program written by Falunas (1996) which 
models the galaxy as a circular, differentially rotating flat 
disk, and assumes no a priori form for the shape of the 
rotation curve. The program fits the velocity data in a 
series of elliptical annuli by adjusting several parameters: 
central position XQ,yo, position angle 0, inclination z, and 
rotation speed Vc within each annulus. Data near the mi- 
nor axis are excluded since they have almost no signal due 
to projection. In practice the kinematic center is not well 
constrained (see also Begeman 1987, Falunas 1996) and 
we fixed the central position to be the photometric center 
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Fig. 7.— VLA map of NGC 4123 and NGC 4116 
The 21 cm surface brightness of NGC 4123 and NGC 4116 is plotted as grayscale, with velocity contours overlaid. North is up and west is to the right. 
The noise increases in the south and west as these regions are at the edge of the VLA primary beam. The surface brightness grayscale runs from Nh — 
to 1.5 X 10 cm~ . The contours are of heliocentric velocity from 1225 to 1425 km s~ at intervals of 25 km s~ , and each galaxy has the receding 
side in the northwest. 



of the galaxy. 

Rotation curve fits to two-dimensional velocity fields 
have degeneracies between the parameters, and are in- 
fluenced by non-circular motions in the disk and warps. 
We examined the variation of </), i, and with radius. 
Inside the bar, the parameters are "twisted" by the non- 
circular motions in the bar. Outside the bar radius, from 
^ 130" to 250", (p, i, and Vc change only slowly and we 
adopted median values of </> and i from this region, yielding 
(p = -33 ± 2° and i = 45 ± 4°. Vc{R) and i{R) are covari- 
ant and show artificial excursions at certain radii, due to 
spiral arms crossing the annuli (see Palunas 1996). How- 
ever, by averaging over a radial range we obtain a reliable 
estimate of the parameters 4> and i. 

We then extract the rotation curve Vc{R) in annuli us- 
ing the adopted values of xq, yo, 0, and We chose to use 
global values of </> and i rather than doing a tilted-ring ex- 
traction, since tilted-ring fits do not work well on galaxies 
of i ~ 45° {e.g. Begeman 1987) and are also thrown off by 
spiral arms. The rotation curve extraction was spatially 



oversampled, in part to avoid further smoothing caused by 
binning into annuli. The residual map in the right panel 
of Figure ^ shows that the model produces a good repre- 
sentation of the data, justifying the use of global and 
i. Outside the bar region, residuals are generally below 10 
km s~^. 

The rotation curve computed from the H I data is shown 
as the points and errorbars in each panel of Figure ^. The 
approaching and receding sides of the galaxy are plotted 
separately. The rotation curves for the two sides differ by 
10-20 km s~^ at large radii; the difference is caused by 
spiral non-circular streaming motions. The right panel of 
Figure |^ shows the residual velocity field after subtraction 
of the axisymmetric rotation curve (using an average of the 
two sides; the residual map is similar if the approaching 
and receding sides are subtracted separately). There are 
strong antisymmetric residuals in the bar region, and a 
spiral-shaped residual from the outer spiral pattern. 

Outside the bar radius of 5.6 kpc, it is reasonable to de- 
rive a rotation curve assuming that the gas is on circular 
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Fig. 8. — Mass models for the rotation curve of NGC 4123 
The contributions of the stellar disk, gas disk, and best-fit halo to the rotation curve of NGC 4123 are shown for models with low and high mass disks, 
and three different dark halo types. 
Points and crrorf;)ars - observed rotation curve. 
Uppermost sofid line - };)cst-fit fotal rotation curve. 
Dotted line - contribution of stars+gas to the rotation curve. 
Lower solid line - contribution of the stellar disk. 
Dashed line — contribution of the best-fit halo. 
Dot-dashed line — contribution of the gas disk. 

The left column has disk M/ L — 1.0, and the right column has disk M/L = 2.25. The top row uses isothermal halos, the middle row uses NFW halos, 
and the bottom row uses Moore-type halos. 



orbits. Inside that radius, the derived rotation curve is af- 
fected by noncircular motions and the error bars increase 
because the data are not fully consistent with circular or- 
bits and because the H I disk has a central hole. 

The maximum of the rotation curve is 145 km s~^, while 
the flat part of the rotation curve is at 14 = 130 km s~^. 
Comparing the rotation width and total / magnitude to 
the TuUy-Fishcr relation of Giovanclli et al. (1997) sug- 
gests that the galaxy is ovcrluininous by 0.2 magnitudes, 
which is consistent with the scatter in that TF relation. 
The difference could also be caused by a 10% overestima- 
tion of the distance to NGC 4123. which would require a 
10% increase to the M/L values we discuss below. 

5.4. The spiral pattern speed 

The spiral velocity residual is a one-armed pattern that 
changes sign several times along the arm. An m = 1 ve- 



locity residual is expected inside corotation of an m = 2 

spiral pattern (Canzian 1993), but an m = 3 pattern is 
expected outside corotation, which suggests that nearly 
all of the H 1 disk of NGC 4123 is inside corotation of 
the spiral. However, the m = 1 spiral residual is not ex- 
actly as in Canzian's model, where the sign of the residual 
does not change along the arm. Alternatively, if the resid- 
ual is interpreted as a streaming motion, it is consistent 
with inward streaming along the arm until the the last 
sign change, at 200" along the major axis southeast of 
the galaxy center. Past this point the positive residual on 
the approaching, far side of the galaxy indicates outward 
streaming. The change from inward to outward streaming 
occurs at the corotation radius, which is then at 21 kpc. 
Again most of the H I disk is inside corotation. 

In Paper 11, we show that the corotation radius of the 
bar is well within the H I disk, at 5 to 9 kpc. Thus the bar 
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and spiral pattern have different corotation radii and dif- 
ferent pattern speeds, even though the spiral arms appear 
to emerge from the ends of the bar. These dynamically 
determined pattern speeds confirm the suggestion of Sell- 
wood & Sparke (1988) that bar and spiral pattern speeds 
in the same galaxy can be significantly different. 

6. MASS MODELS 
6.1. Fitting the rotation curve 

We fit axisymmetric mass models made up of the stel- 
lar disk, gas disk, and a dark halo to the rotation curve. 
Because they are axisymmetric, these models neglect the 
bar distortions both in the light and in the velocity field, 
which are strong only inside the bar radius of 5.6 kpc. 
Our goal is to determine the halo required to fit the outer 
rotation curve for a range of values for the disk mass-to- 
light ratio. We also want to determine the maximum disk 
mass-to-light ratio, the highest disk M/L consistent with 
the rotation curve data and a non-hollow dark halo. 

We calculate the rotation curve of the axisymmetrized 
stellar disk from the /-band image, assuming a uniform 
M/L. Our procedure is described in more detail in Pa- 
per II. Briefly, we remove foreground stars and rectify the 
image to face-on, using the position angle and inclination 
obtained from the H I observations. We then average the 
image in circular annuli about the center to produce an 
axisymmetrized surface brightness distribution. We as- 
sume that the stellar disk has a vertical distribution of 
the common form p{z) oc sech^(z/zo) and a scale height of 
zq = 400 pc, similar to that of the disk of the Milky Way 
(Mihalas & Binney 1981). Following Quillcn et al. (1994), 
we calculate the gravitational accelerations as a function 
of position, summing the contributions of each pixel ef- 
ficiently using FFT techniques (Hockney 1965). The ac- 
celeration as a function of radius determines the stellar 
contribution Vdisk{R) to the rotation curve Vc{R). 

The rotation curve due to the stellar disk for an assumed 
AI/L is plotted as the lower solid line in the panels of Fig- 
ure H It is plotted out to the end of the optical image 
only, though its contribution is accounted for outside this 
radius. Its slightly unusual shape, rising quickly, leveling 
off, then rising again, is due to the shoulder in the ax- 
isymmetrized surface brightness profile (Figure |^) in the 
bar region. 

We calculate the rotation curve due to the gas disk in a 
similar fashion. We assume that the H I disk is thin and 
flat and rectify it to face-on, average the column density 
azimuthally, and scale the mass up by 1.4 to correct for 
helium. The atomic gas contribution to the rotation curve, 
shown as the dot-dashed line in Figure B, is everywhere a 
rather small fraction of the observed velocity, and inside 
7 kpc it makes a very small negative contribution to Vc^ 
(i.e. outward centripetal acceleration) due to the central 
hole in the atomic gas density. The dotted line in Figure || 
is the combined contribution of the stellar and gas disks - 
the contribution of visible, non-dark matter to the rotation 
curve. 

We do not have CO imaging data for NGC 4123, and 
thus cannot make a detailed calculation of the contribu- 
tion of its molecular gas to the rotation curve. Single-dish 
CO observations with the SEST and IRAM 30-m imply a 
molecular gas mass of 4.0 X 10^ M0 within a 23" beam, and 



5.1 X 10® Mq within a 45" beam (Chini, Kruegel & Steppe 
1992). The contribution to the rotation curve from the 
molecular gas is then 37 km s^^ at 1.25 kpc and 30 km 
at 2.5 kpc. The molecular gas may at most fill in the hole 
in the atomic gas distribution, but its contribution to the 
overall rotation curve is negligible. 

The dark halo needed to fit the observed rotation curve 
is specified by two parameters: a characteristic density and 
radius. We performed least-squares fits for these parame- 
ters over a range of disk M / L and three different models 
for the dark halo: a pseudo-isothermal model with a con- 
stant density core and r~^ envelope: 
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and two broken power-law models, the "NFW profile" and 
"Moore profile" found in A^-body simulations of dark mat- 
ter halo formation (Navarro, Frenk, & White 1996; Syer 
& White 1998; Moore et al. 1999). The NFW-type model 



has an inner r 
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and outer r ^ dependence, and break 
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This pseudo-NFW halo model is not precisely an "NFW 
model," because the NFW family of halos is a one-para- 
meter family: Navarro, Frenk & White (1996) found a 
relation between the characteristic density and radius for 
halos formed in any given cosmology, whereas we allow 
these parameters to vary independently in order to fit the 
rotation curve. An NFW-type profile can equivalently be 
uniquely expressed as a function of concentration index c 
and virial velocity V200 (Navarro et al. 1996). The Moore- 
type profile is similar to the NFW halo but has an inner 
cusp slope of 



-1.5. 
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Halos with this profile are found in high-resolution A^-body 
dark matter simulations {e.g. Moore et al. 1999). 

These expressions specify spherical mass distributions, 
although the rotation curve constrains only the centripetal 
acceleration in the disk plane; for any degree of halo flat- 
tening, a volume density distribution could be devised to 
yield the same rotation curve. The functional forms used 
are analytically convenient, but are not intended to favor 
specific models for halo formation; we have also not at- 
tempted to modify the halo profiles for compression by 
disk collapse, a point we discuss further in Paper II. 

The fits are not very sensitive to the functional form 
adopted: all three types of halo model fit the data well 
for disk M/L < 2.25, although for high-mass disks with 
M/L > 2.0, the cuspy NFW and Moore halos are forced 
to very low scale densities and large scale radii in order to 
contribute to the rotation curve at large radius but avoid 
exceeding the observations at i? ~ lOkpc. 

We constructed models over the range M/L = 1.0 to 
3.0 with the three halo profiles. The models and best-fit 
halo parameters are listed in Table ^. For all profile types, 
the fitted halo density and radius are strongly covariant 
in a way which keeps the halo contribution to the outer 
rotation curve roughly constant, although the actual error 
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ellipse on the fitted parameters is quite narrow - an exam- 
ple is shown in Paper II. We use the halos fitted in these 
mass models in the fluid dynamical simulations described 
in Paper II. 

Columns 7 and 8 of Table |^ show the results of fitting 
NFW-type profiles by fitting for c and ¥200- The best- fit 
halo profile p{r) is the same as that obtained by fitting 
for pg and , but c and V200 &re much less covariant than 
Ps and Vs and are tightly constrained for a given M/L. 
As disk M/L is increased, the halo contribution inside the 
optical radius decreases greatly and the radius at which 
the halo rotation curve peaks moves out, requiring a lower 
concentration index c. In order to maintain the necessary 
halo contribution at large i?, the virial velocity V200 must 
increase as c decreases; the virial radius r2oo increases with 
V200- For M/L > 2.0, the halo rotation curve peaks be- 
yond the end of the H I data at 47 kpc. 

6.2. Discussion 

We plot rotation curve decompositions for AI/L of 1.0 
and 2.25 in Figure |. For M/L < 2.25, aU mass models fit 
the observed rotation curve well and the small differences 
in the values of reduced x^, given in Table § are not sig- 
nificant, given that the true shape of the halo radial profile 
is unknown and could be adjusted slightly to improve the 
fit. 

Models with too high a disk M/L cannot fit well, how- 
ever, because the fitted rotation curve exceeds the ob- 
served curve. The stellar disk alone accounts for the entire 
rotation curve at 8 kpc when M/L — 2.5 and the model is 
an acceptable fit since it only marginally exceeds the data 
when the isothermal halo is included. For higher M/L, the 
total rotation curve is too high, and the fit is unacceptable. 
For an isothermal halo, therefore, the maximum disk M/L 
is about 2.5. NFW and Moore-type halos, which are more 
centrally concentrated, contribute more to the peak of the 
rotation curve, and for them the maximum disk M/L is 
slightly less, about 2.25. For maximum disk models, the 
NFW and Moore halos have low densities and large scale 
radii. 

6.3. Radial variation in M/L 

It is conventional to assume a constant M/L when fit- 
ting rotation curves, but it is not guaranteed, although 
there is some evidence for constant M/L (van der Kruit 
1981). There is a color gradient in NGC 4123, as discussed 
in Section ^.3| , which could be caused by a change in the 
disk stellar population. To test the effect of a radial vari- 
ation in M/L, we compared the stellar disk contribution 
to the rotation curve for a constant M/L model and for 
models in which the Al / L declines outside the bar radius. 
The stellar population is likely to be radially mixed inside 
the bar radius, and in Paper II we will determine the stel- 
lar M/L inside the bar radius, so the important question 
is whether the rotation curve decomposition is strongly af- 
fected by a disk M/L which differs inside and outside the 
bar radius. Figure O shows the disk contribution for three 
models: a constant M/L — 2.25 model, and two models 
with a steep decline in M/L outside the bar radius. 

Even these models with a rather exaggerated decline in 
M/L only have a small effect on the disk contribution to 
the rotation curve. Reasonable variations in M/L there- 
fore do not substantially increase the amount of dark mat- 
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Fig. 9. — Effect of varying M/L on disk rotation curve 
The disk contribution to the rotation curve is shown for three models. 
Solid line: M/L— 2.25, constant; dashed line: M/L dechning outside 6 
kpc, to 0.5 at 20 kpc; dotted line: M/L declining outside 5 kpc, to at 
20 kpc. 



ter required by the rotation curve. An increase in disk 
M/L with radius, as could happen if the outer disk is 
bluer because it is old and metal-weak, would of course 
reduce the required dark matter, but again the effect will 
be small. 

The effect of varying M/L is small because at the bar 
radius of 5.6 kpc, the disk contribution to the rotation 
curve is already nearing its peak. The decline from the 
peak of a rotation curve is relatively gradual, whether or 
not the mass distribution is cut off. A declining M/L out- 
side the bar radius cannot make the rotation curve decline 
fast enough to change the inferred halo parameters signif- 
icantly - in fact, a declining M/L makes the disk part of 
the rotation curve peak up by a tiny amount, because it 
behaves like a disk truncation. 

7. CONCLUSIONS 

NGC 4123 is a strongly barred galaxy with a rather boxy 
bar containing two components, both of which are more 
rectangular than ellipsoidal. The bar has the straight off- 
set dust lanes characteristic of strong bars (Athanassoula 
1992b), although they are less obvious in the / image than 
in bluer bandpasses. The outer disk of the galaxy shows 
a roughly exponential decline of surface brightness with 
radius to the limit of the data, at / ~ 28; B,V ^ 29. 
The disk extends to at least 10 exponential disk scale 
lengths; there is no sign that it is truncated, in agreement 
with 1 mag deeper measurements of NGC 5383 (Barton 
& Thompson 1997) but in sharp contrast to the conven- 
tional result that disks are truncated at 3-5 scale lengths 
(cf. van der Kruit & Searle 1981a,b, 1982; van der Kruit 
1989; Barteldrees & Dettmar 1994). Blue knots in the 
outer spiral arms indicate that stars are forming well out- 
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0.72 
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81.0 ±0.8 


0.522 


19.0 
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1.75 


36.1 


2.24 


0.84 


0.92 


14.2 


0.77 


6.76 ±0.06 


82.6 ±0.9 


0.237 


26.9 


0.83 


2.00 


12.7 


3.76 


0.87 


0.39 


21.4 


0.85 


4.35 ±0.04 


87.5 ± 1.0 


0.0802 


45.4 


0.96 


2.25 


4.68 


6.33 


0.93 


0.13 


38.2 


1.02 


2.09 ±0.03 


93.0 ± 1.2 


0.0399 


61.0 


1.33 


2.50 


1.91 


10.4 


1.11 


0.043 


69.1 


1.40 


1.30 ±0.01 


116.3 ±1.7 


0.0101 


128.3 


1.93 


2.75 


0.88 


16.9 


1.63 


0.011 


185.7 


2.10 


0.645 ± 0.005 


142.8 ±2.2 


0.0063 


156.3 


3.09 


3.00 


0.46 


27.1 


2.75 


0.0038 


381.7 


3.46 


0.255 ± 0.002 


151.1 ±2.9 


0.0040 


179.4 


4.77 



Table 4 

Mass models: best-fit halo parameters 



Note. — The best-fit halo parameters for disk mass-to-light ratios 
from 1.0 to 3.0. Column 1 is the disk M/L. Columns 2, 3, and 4 are 
for isothermal halos; columns 5, 6, and 7 are for NFW-type power-law 
halos; and columns 10, 11, and 12 are for "Moore-type" power-law 
halos. The densities are in 10~^ Mq pc~^, and the radii are in kpc. 
Columns 8 and 9 are the parameters and la errors for the best fit 
NFW-type halos when the fitted parameters are the concentration 
index c and virial velocity V200 ■ 



side i?25- 

Fabry-Perot observations of the velocity field show strong 
non-axisymmetric motions inside the bar radius, and large 
velocity jumps on the leading edges of the bar, at the lo- 
cation of the dust lanes. These support the interpretation 
of the dust lanes as the location of shocks in the gas flow. 
The velocity jumps are offset and several velocity contours 
make a sharp jog as they passes through the center of the 
galaxy, suggesting the presence of an inner Lindblad res- 
onance. The strength of the non-axisymmetric motions 
decreases rapidly beyond the bar radius. 

21 cm observations of NGC 4123 show that it has an ex- 
tended gas disk in circular rotation outside the bar radius, 
with spiral structure visible in the intensity and weakly 
in the kinematics. The kinematics indicate that corota- 
tion of the spiral pattern is at ~ 4 times the bar radius; 
when combined with the result from Paper II that the bar 
is fast-rotating, this implies that the bar and spiral have 
very different pattern speeds. There is no evidence for 
a warp. The rotation curve derived from these observa- 
tions extends to 4 optical radii (4i?25) and requires the 
presence of a dark matter halo in addition to the lumi- 
nous disk. The maximum height of the rotation curve is 
~ 145 km s"^; the height of the flat part of the rotation 
curve is 130 km s~^. 

Mass models that include a stellar disk and dark halo to 
the H I rotation curve yield acceptable fits, in a sense, 
for /-band disk mass-to-light ratios of 2.5 or less. Model 
halos with either isothermal or power-law density distribu- 
tions fit the data well. The axisymmetric rotation curve 
does not contain enough information to constrain the disk- 
halo decomposition further. Both heavy and light disk 
models, and both isothermal and power-law dark halos, 
yield acceptable fits to the rotation curve. 

However, there is extra information in the non-circular 
gas streaming motions seen in the Fabry-Perot data. In 
Paper II, we use this information to show that the disk of 
NGC 4123 is close to maximum, with M/Lj = 2.25 the 



preferred value. 

We thank the staff at CTIO, the VLA, and Las Cam- 
panas for their excellent support during these observa- 
tions. BJW thanks Karl Gebhardt and Povilas Pahmas 
for Fabry-Perot reduction software and helpful comments. 
We thank Jason Prochaska for a reading of the manuscript 
and the anonymous referee for helpful comments. This re- 
search was supported in part by NSF grant AST 96/17088 
and NASA LTSA grant NAG 5-6037 to JAS, by NSF grant 
AST 96-19510 to TEW, and NSF grant AST 96-17177 to 
JvG. BJW has been supported by a Barbara McClintock 
postdoctoral fellowship from the Carnegie Institution of 
Washington. Operation of the RFP is partially supported 
by CTIO. 



REFERENCES 



Athanassoula, E. 1992a, MNRAS, 259, 328 
Athanassoula, E. 1992b, MNRAS, 259, 345 

Athanassoula, E., Bosnia, A., & Papaioannou, S. 1987, A&A, 179, 23 
Babcock, H.W. 1939, Lick Observatory Bulletin, 19, 41 
Bahcall, J.N., & Casertano, S. 1985, ApJL, 293, L7 
Barteldrees, A. & Dettmar, R.-J. 1994. A&AS. 103. 475 
Barton, I.J., & Thompson, L.A. 1997, A.I. 114. 655 

Beck, R., Ehle, M., Shoutenkov, V., Shukurov, A., & Sokoloff, D. 1999, 

Nature, 397, 324 
Beers, T.C., Flynn, K., fc Gebhardt, K. 1990, A,T, 100, 32 
Begeman, K. 1987, PhD thesis, Groningen 
Bosma, A. 1981, AJ 86, 1825 

Bosnia, A. 1999, in Galaxy Dynamics, ed. Merritt, D.R., Valluri, M., &: 

Sellwood, J. A. (San Francisco: A.S.P.), 339 
Bottema, R. 1993, A&A, 275, 16 
Bottema, R. 1997, A&A, 328, 517 

Briggs, D.S. 1995, PhD thesis. New Mexico Institute of Mining and 

Technology 
Briggs, F.H. 1990, ApJ, 352, 15 
Bruzual A., G. & Chariot, S. 1993, ApJ, 405, 538 
Burstein, D., & Heiles, C. 1982, AJ, 87, 1165 
Byun, Y.I., Freeman, K.C. & Kylafis, N.D. 1994, ApJ, 432, 114 
Cardelli, J. A., Clayton, G.C. & Mathis. J.S. 1989, ApJ, 345, 245 
Chariot, S., Worthey, G., & Bressan, A. 1996, ApJ, 457, 625 
Chini, R., Kruegel, E. & Steppe, H. 1992, A&A, 255, 87 
Courteau, S., & Rix, H.-W. 1999, ApJ, 513, 561 
de Jong, R. 1996a, A&A, 313, 45 



16 



Disk and Halo Mass of NGC 4123. I. 



do Jong, R. 1996b, A&A, 313, 377 

dc Vaucoulcurs, G. , et a/. 1991, Third Reference Catalog of Bright 

Galaxies (RC3), (New York: Springer- Verlag) 
Debattista, V.P. & Sellwood, J. A. 1998, ApJL, 493, L5 
Duval, M.F., & Athanassoula, E. 1983, A&A, 121, 297 
Efstathiou, G., Lake, G. & Negroponte, J. 1982, MNRAS, 199, 1069 
England, M.N. 1989, ApJ, 337, 191 

Ferguson, A.M.N., Wyse, R.F.G., Gallaglier, J.S. & Hunter, D.A. 1998, 

ApJ, 506, L19 
Freeman, K.C. 1970, ApJ, 160, 811 

Freeman, K.C. 1992, in Physics of Nearby Galaxies: Nature or Nurture?, 

od. Thuan, T.X., Balkowski, C, & Van, J.T.T. (Gif-sur-Yvotte: 

Editions Frontieres), 201 
Giovanelli, R., Haynes, M., Salzer, J. Wegner, G.. da Costa. L., & 

Freudling, W. 1994, AJ, 107, 2036 
Giovanelli, R., Haynes, M.P., Herter, T., Vogt, N.P., Wegner, G., Salzer, 

J. J., Da Costa, L.N., & Freudling, W. 1997, AJ, 113, 53 
Graham, J. A. 1982, PASP, 94, 244 

Hockney, R.W. 1965, J. Assoc. Comp. Mach., 12, 95 

Huchtmeier, W.K., & Richter, O.-G. 1989, A General Catalog of H 1 

Observations of Galaxies, (New York: Springer) 
Kalnajs, A. 1983, in Internal Kinematics and Dynamics of Disk Galaxies, 

lAU Symposium 100, ed. Athanassoula, E. (Dordrecht: Reidel), 87 
Kennicutt, R. C. , Jr. 1989, ApJ, 344, 685 
Landolt, A. U. 1992, AJ, 104, 340 

Lindblad, P.O., Hjelm, M., Hogbom, J., Jorsater, S., Lindblad, P.A.B., 

Santos-Lleo, M. 1996a, A&AS, 120, 403 
Lindblad, P.A.B., Lindblad, P.O., & Athanassoula, E. 1996b, A&A, 313, 

65 

Long, K. 1991, PhD thesis, Princeton University 

MacAlpine, G.M., & WilUams, G.A. 1981, ApJS, 45, 113 

Mihalas, D., & Binney, J. 1981, Galactic Astronomy, (New York: Freeman) 

Moore, B., Quinn, T., Governato, F., Stadel, J. & Lake, G. 1999, 

MNRAS, 310, 1147 
Moore, E.M., & Gottesman, S.T. 1995, ApJ, 447, 159 
Morrison, H.L., Boroson, T.A., & Harding. P. 1994, AJ, 108, 1191 
Navarro, J.F., Frcnk, C.S., & White, S.D.M. 1996, ApJ, 462, 563 
Navarro, J.F. 1998, preprint Eistro-ph/980708d 
Ondrechen, M.P., & van derBulst, J.M. 1989, ApJ, 342, 29 
Ondrechen, M.P., van der Hulst, J.M., & Hummel, E. 1989, ApJ, 342, 39 
Opik, E. 1922, ApJ, 55, 406 

Palunas, P. 1996, PhD thesis, Rutgers University 

Palunas, P., & Williams, T.B. 2000, ApJ, submitted 

Pence, W.D., & Blackman, CP. 1984a, MNRAS, 207, 9 

Pence, W.D., & Blackman, CP. 1984b, MNRAS, 210, 547 

Finer, B.C., Stone, J.M., & Teuben, P.J. 1995, ApJ, 449, 508 

Prendcrgast, K.H. 1966, in Theory of Orbits in the Solar System and 

in Stellar Systems, lAU Symposium 25, cd. Contopoulos, G., (New 

York: Academic), 60 
Prendcrgast, K.H. 1983, in Internal Kinematics and Dynamics of Galaxies, 

lAU Symposium 100, ed. E. Athanassoula, (Dordrecht: Reidel), 215 
Press, W.H., Flannery, B.P., Teukolsky, S.A. & Vetterling, W.T. 1992, 

Numerical Recipes, (Cambridge: Cambridge U.P.) 
Quillen, A.C., Frogel, J. A., & Gonzalez, R.A. 1 994, ApJ, 437, 162 
Quillen, A.C., Frogel, J. A., Kuchinski, L.E. & Tcrndrup, D.M. 1995, AJ, 

110, 156 

Regan, M.W., Vogel, S.N., & Teuben, P.J. 1995, ApJ, 449, 576 
Regan, M.W., Vogel, S.N., & Teuben, P.J. 1997. ApJL, 482, L143 
Rhoads, J.E. 1998, AJ, 115, 472 
Roberts, M.S. 1969, AJ, 74, 859 

Rubin, V.C, Ford, W.K., Jr., & Thonnard, N. 1978, ApJ, 225, L107 
Sandage, A. 1961. The Hubble Atlas of Galaxies, (Washington: Carnegie 

Inst, of Washington) 
Sandage, A., & Bedke, J. 1988, Atlas of Galaxies: Useful for Measuring 

the Cosmological Distance Scale, (Washington: NASA, US GPO) 
Schommer, R.A., Bothun, G.D., WiUiams, T.B., & Mould, J.R. 1993, 

AJ, 105, 97. 
Schwarz, M.P. 1984, MNRAS, 209, 93 
Sellwood, J. A. & Balbus, S.A. 1999, ApJ, 511, 660 
Sellwood, J. A. & Sparke, L.S. 1988, MNRAS, 231, 25P 
Syer, D., Mao, S., & Mo, H.J. 1999, MNRAS, 305, 357 
Syer, D. & White, S.D.M. 1998, MNRAS, 293, 337 

Taylor, C.L., Brinks, E., Grashuis, R.M., & Skillman, E.D. 1995, ApJS, 
99, 427. 

Tully, R.B., & Fisher, J.R. 1977, A&A, 54, 661 

van Albada, T.S., Bahcall, J.N., Begeman, K., & Sancisi, R. 1985, ApJ, 
295, 305 

van Albada, T.S., & Sancisi, R. 1986, Phil. Trans. Roy. Soc. Lond. A, 
320, 447 

van Albada, T.S., & Sanders, R.H. 1982, MNRAS, 201, 303. 
van der Kruit, P.C. 1981, A&A, 99, 298 
van der Kruit, P.C, & Searle, L. 1981a, A&A, 95, 105 
van der Kruit, P.C, & Searle, L. 1981b, A&A, 95, 116 
van der Kruit, P.C, & Searle, L. 1982, A&A, 110, 61 

van der Kruit, P.C. 1995, in lAU Symposium 164, Stellar Populations, 
ed. van der Kruit. P.C. & Gilmore, C, (Dordrecht: Kluwcr), 205 

van der Kruit, P.C. 1989, in The World of Galaxies, cd. Corwin, H.G. & 
Bottinelli. L. (New York: Springer), 256 

Wciner, B.J., & Sellwood, J. A. 1999, ApJ, 524, 112 

Weiner, B.J., Sellwood. J. A., & WilUams, T.B. 2000, ApJ in press (Paper 
II) 

Worthey, G. 1994, ApJS, 95, 107 



